DISSIPATIVE BOUSSINESQ SYSTEM OF EQUATIONS 
IN THE BENARD-MARANGONI PHENOMENON 



R. A. Kraenkel, S. M. Kurcbart, J. G. Pereira 
Instituto de Ftsica Tedrica 
Universidade Estadual Paulista 
■ Rua Pamplona 145, 

o\ ' 

00 
(N 



O 

C/3 



Ph. 



01405-900 Sao Paulo SP - Brazil 

M. A. Manna 
Laboratoire de Physique Mathematique 



^ ' Universite de Montpellier II 

(N 



34095 Montpellier Cedex 05 - France 



o 
o 

O ■ (February 9, 2008) 

cn 



Abstract 



By using the long- wave approximation, a system of coupled evolution 



^ ' equations for the bulk velocity and the surface perturbations of a Benard- 



Marangoni system is obtained. It includes nonlinearity, dispersion and dissi- 
pation, and it can be interpreted as a dissipative generalization of the usual 
Boussinesq system of equations. As a particular strictly dissipative 

version of the Boussinesq system is obtained. Finally, some speculations are 
made on the nature of the physical phenomena described by this system of 
equations. 

PACS: 47.20.Bp ; 47.35.+i 
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I. INTRODUCTION 



Several recent works |]T]-§1 have dealt with the study of oscillatory instabilities in sys- 
tems of the Rayleigh-Benard and Benard-Marangoni type. Considering a shallow fluid layer 
bounded below by a plane stress-free perfect conducting plate, and with a deformable upper 
free surface where a constant heat flux is imposed, the linear stability analysis indicates that 
oscillatory motion sets in when a certain critical combination of the Rayleigh and Marangoni 
numbers is attained Subsequently, it was shown that, when the surface-tension can be 
disregarded and the Rayleigh number of the system is at i? = 30, appropriate surface deflec- 
tions governed by the Korteweg-de Vries (KdV) equation appear . For a two-dimensional 
surface, the same deflections turn out to be governed by the Kadomtsev-Petviashvili (KP) 
equation 0. Thereafter, these results have been extended to the Benard-Marangoni sys- 
tem 0,^ , where buoyancy is discarded. In this case, appropriate surface disturbances prop- 
agate according to the KdV equation when the Marangoni number of the system is at 
M = —12. Further extentions to a double-diffusive systems [Q, and to a temperature 
depending viscosity have also been established. The physical mechanism behind such 
sustained excitations is the balance, at the critical point, between the energy dissipated by 
viscous forces and that released either by buoyancy, or by a temperature depending surface 
tension. In both cases, out of the critical points, appropriate surface excitations are governed 
by the Burgers equation [^,||. 

The KdV equation is well known to govern long surface-waves in inviscid shallow flu- 
ids It corresponds to a situation where nonlinearity and dispersion compensate each 
other, making possible the existence of coherent wave-structures, like the solitary-wave. 



The reductive perturbation method of Taniuti ||I0[, based on the concept of stretching, and 
in which waves propagating in only one direction are sought from the beginning, is a com- 
mon approach to study long waves in shallow water. From the various possible stretchings 
of the coordinates, different evolution equations (linear wave, breaking wave, KdV) may 
emerge as governing surface disturbances. However, there is an alternative approach to the 
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theory of long waves in shallow water, which is based on perturbative expansions in two 
small parameters [^. One of them measures the smallness of the amplitude perturbation, 
and the other is a measure of the longness of the wave-length perturbation. This approach, 
as an intermediate step, and from different possible relations between the two perturba- 
tive parameters, yields different systems of evolution equations (nonlinear shallow water, 
Boussinesq) describing superimposed waves propagating in opposite directions. Only when 
specializing to waves moving in a given direction, equations like breaking wave and KdV 
show up. A perturbative scheme of this kind has not been used to study surface excitations 
in Benard systems. Consequently, for these systems, the more general evolution equations, 
wherefrom Burgers and KdV equations are obtained, have not been found. 

In this paper, instead of using the reductive perturbation method of Taniuti [0 , we will 
proceed through a perturbation scheme for the Benard- Mar angoni system based on two per- 
turbative parameters, leading to the so called long-waves in shallow water approximation. 
By this way, a new system of coupled evolution equations will be found, involving the fluid 
velocity and the free-surface displacement. This system may be interpreted as a dissipative 
generalization of the Boussinesq equations. When the Marangoni number assumes the criti- 
cal value M = —12, and a certain relation between the perturbative parameters is assumed, 
it reduces, as we are going to see, to the usual Boussinesq equations. A further restriction 
to waves moving in only one direction leads to the KdV equation. On the other hand, out 
of the critical point, and assuming a different relation between the perturbative parameters, 
the dissipative generalization of the Boussinesq equations reduce to a strictly dissipative 
version of the Boussinesq equations. In this case, a restriction to waves moving in only one 
direction will lead to the Burgers equation. Although we have restricted ourselves to the 
Benard- Marangoni system, an extension to the Rayleigh-Benard case, where buoyancy is 



predominant, may also be established with similar results |Tl]]. Finally, through a linear 
analysis, a speculation on the nature of the physical solutions to the strictly dissipative 
Boussinesq system of equations is made. 
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II. BASIC EQUATIONS AND BOUNDARY CONDITIONS 



We consider a fluid bounded below by a plane, stress-free, perfect termally conducting 
plate at 2; = 0, and above by a deformable one- dimensional free surface which, at rest, lies at 
z = d. The depth d will be supposed to be small enough so that buoyancy can be neglected 
when compared to the effects coming from the surface tension dependence on temperature. 
In other words, we will be dealing with a Benard-Marangoni system. The equations that 
describe such a system are 

V.v = 0, (1) 

dv ^ ^ 

= -Vp + /^V\+(7P, (2) 

f^^v^r, (3) 

where ^ = ^ + "^-V is the convective derivative, v = (u,0,w) is the fluid velocity, and p 
is the pressure. The density p, the viscosity fi, and the thermal diffusivity coefficient k are 
supposed to be constant. The surface tension r will be assumed to depend linearly on the 
temperature: 

r = ro[l-7(T-To)], (4) 

where 7 is a constant, and tq, Tq are reference values for the surface tension and the tem- 
perature, respectively. 



The boundary conditions on the upper free surface z = d + ri{x, t) are |jT2| 



r]t + ur]^ = w , (5) 



{P -Pa)-J^ {i^z + U^VI - VxUz - VxWcc) = -^Vxx , (6) 
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r],T, -T, = jN, (8) 

where F is the normal heat flux, k is the thermal conductivity, pa is a pressure exerted on 
the upper surface, all of them supposed to be constant, and = (1 + ril) '^. Subscripts 
denote partial derivatives with respect to the corresponding coordinate. 

On the lower plane z = 0, we assume that the sliding resistance between two portions 
of the fluid is much greater than between the fluid and the plane |T3[ , implying a stress- free 
lower surface: 

w = u^ = 0. (9) 
Moreover, we will assume that the lower plate is at a constant temperature T = Th. 



III. PERTURBATIVE SOLUTION: EVOLUTION EQUATIONS 

The static solution to the above equations is given by 

Ps=Pa- pgiz - d) , 

F 

Ts = To- —{z - d) . 

We consider now perturbations from this quiescent state. The horizontal and vertical length 
scales of these perturbations are supposed to be I and a, respectively. Then, we define two 
small parameters 

which will be used to order the expansions. Before proceeding further, however, it is con- 
venient to write the equations, boundary conditions and static solutions in a dimensionless 
form. This is done by taking the original variables (primed) to be 

x' = Ix z' = dz t' = —t (lla) 

Co 

rj = ar] u= — u w = — -w, (lib) 
Co Cod 



where Cq = gd. Furthermore, four dimensionless parameters appear: the Pradtl number 
a — ii/pn; the Reinolds number R — codp/ii; the Bond number B — pgd^/ro, and the 
Marangoni number M = ■yPd^To/kKp. 

To obtain the nonhnear evolution of the surface perturbations in the shallow water theory, 
we expand all variables in powers of z, keeping the terms that will contribute to the evolution 
equations up to orders e and 5^. Despite laborious, this procedure is straightforward, and 
for this reason we will omit the details, specifying in the text only the general guidelines, 
and giving a brief summary of the results in the Appendix. To start with, we make the 
expansions 



oo oo 



ra=0 n=0 

where and w„ are both functions of x and t. Then, substituting them in Eq.(l), we get 
the relation 



r-2 '^^'nx 



n + l 

Using the boundary conditions at z = 0, it is easy to show that 

Wo = ty2 = ■w^4 = ■■■ = , 
U\ = = = ... = Q . 

Then, using the expansion 



n=0 

in Eq. (2), it is possible to obtain it2, ti4, Uq and p2 in terms of Uq and po only. The other 
components of the expansions will contribute to orders higher than e and 5^, and therefore 
they can be neglected. 

Next, expanding the temperature according to 

oo 
n=0 

and using Eq. (3) with the corresponding boundary conditions, we can see that 

9r\ — 9-y — 9a — Of, — ... — Q . 



Moreover, we can also obtain expressions for ^3 and $5 in terms of uq and po only. Now, 
Eq.(6) yields po in terms of Uq and rj. Consequently, it is possible to rewrite the w's, p's and 
6''s in terms of uq and r] only. Finally, using the above results in Eqs. (5) and (7), we obtain, 
up to order e and 5^, a coupled system of evolution equations for uq and 77: 



(4cr - 1) 

6 30 ^ ! 



/ 1 ^ \ 

uoxxt -d ' ;g + 30 ■'" ' ^^^^ ^ ^ ' 



(12) 



6^R^ 



Vt + Uox + € (uorj)^ + — {Uoxt + rjxx) - -jUq^xx + {uoxtt + Vxxt) = , 



where 



1 - 



M 



The velocity uq is only the first term in the expansion of u, which is : 

u = Uo + d—^ I Uot + r]x- -^uqxx + —— [um + "nxt) + C [ed, 6 j . 



2 V ^ 

The value averaged over the depth is 

u = U(i + d— {uot + rix)- — 



The inverse is 



R 5"^ 
{rix + ut) + — 



UOxx - T77 i^Ott + Vxt) 

dU 



ioU 



+ O (e5, S^) . 



+ O [eS, S^) 



(13) 



(14) 



Substituting this in Eqs. (12) and (13), and using the lowest order equations into the Uxxt 
term, we obtain (omitting the tilde): 

ut + c r}x + euux ~ ( ^ + ) ""^^ + " ^^^^^ ^ ^ *^1^^) 
?7t + lio; + e(M^)x = , (15b) 



where 



Due to the presence of the Uxx term in Eqs. (15), this system can be considered as a 
dissipative generahzation of the Boussinesq equations. When the Marangoni number is at 
the critical value M = — 12, these equations coincide with the usual Boussinesq system of 
equations 

ut + c^Vx + (^uux + S'^^Vxtt = (17a) 
r]t + Ux + e{ur])x = 0. (17b) 

In this case, by assuming that 6"^ ~ e, and by speciahzing to waves moving, say, to the right, 
we can obtain a relation between u and t]: 

u = CT] - e^T]"^ - 6^^r]xt, (18) 

which is a kind of Riemann invariant [^. Using this in Eq. (15) yields 

3c 2^ 

Vt + crjx + e—T]T]x + d —r]xxx = , (19) 

where now 

Equation(19) is the KdV equation. 

Let us now consider the case M ^ —12. Assuming that 5 ~ e, and neglecting terms of 
order 6^ ~ e^, Eq. (15) becomes: 



ut + c 7]^ + euux ~ ^\^ + = (21a) 

'nt + Ux + e{ur])x = . (21b) 

This is a strictly dissipative version of the Boussinesq system, since the dispersive term is 
not present now. By specializing again to waves moving to the right, we obtain the following 
relation between u and rj: 

u = cr/--r-— 4 + — U,. (22) 



4 ' 2R\ 3 J 



Substituting this relation in Eq. (21), we obtain 



Vt + CT]^ + ^^mx - 77^ 4 + — r/^.^. = , (23) 



2 " 2R\ 3 

which is the Burgers equation. 



IV. THE STRICTLY DISSIPATIVE BOUSSINESQ EQUATION 

The usual Boussinesq system of equations, Eqs. (17), may be transformed into the 



classical Boussinesq equations, also known as dispersive long- wave equations |Tj]. These 



equations have already been shown to be integrable [0. The solutions to KdV and Burgers 



equations have also been extensively discussed in the literature [0. However, the strictly 
dissipative Boussinesq system of equations, Eqs. (21), seems not to have been handled. The 
existence of analytical solutions, therefore, is still an open issue. 

The purpose of this section is to make some speculations on the nature of the physical 
phenomena governed by the strictly dissipative Boussinesq system of equations. To start 
with, we first rewrite Eqs. (21) in a dimensional form: 

Ut + c^gtlx + — adcQUxx = (24a) 
rjt + dux + u^rj + urjr, = , (24b) 



with 



1 ^ M\ 



and given by Eq. (14). For M = —12, and changing the origin of the coordinates 
according to 

ri{x, t) = h{x, t) — d , 



we obtain 



Ut + (?ghx + = (25a) 
ht + Uxh + uhj. = (25b) 
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which, when = 1, coincides with the classical one- dimensional shallow water system of 
equations ||^. A general analytical solution to this system, which might present shocks, has 



already been reported in the literature [|T^]. In fact, by restricting to waves moving to the 
right we obtain the breaking-wave equation [R] 



hf — -\fqdhr + —\f^,hh^ 
2 ^ 2\l d 



0. 



(26) 



On the other hand, for M ^ —12, it is possible to gain some insight on the nature of 
the solutions to Eqs. (24) by making a linear analysis. Neglecting the nonlinear term, Eqs. 
(24) reads: 



ut + (?gri^ - adcQUxx = 
rit + dux = 0. 



(27a) 
(27b) 



Supposing solutions of the form exp{Kx + ut), we obtain the following dispersion relation: 



adc^H? adc^H? 

oj = i ± 

2 2 

There are three cases to be considered: 
(i)When Ac^/a^d^K^ < 1, then 

2c 



-1 + 



4c' 



,2 1 1/2 



(28) 



K > 



ad 



or 



K < - 



2c 
ad 



and the frequency is a pure imaginary number. In this case the system will present either, 
a strictly dissipative or antidissipative behavior, depending on whether the imaginary part 
of uj is negative or positive. 
(ii)When 4c^/a^(i^K^ > 1, then 

2c 2c 
ad ad 

and the frequency is a complex number. In this case the system will present oscilla- 
tory motion combined with either dissipation or anti-dissipation, depending on whether the 
imaginary part of u is negative or positive. For n within this interval, the nonlinear dissi- 
pative Boussinesq system, Eqs. (24), may possibly present a constant-profile opposite two 
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travelling-wave solutions if a compensation mechanism between dissipation and nonlinearity 
takes place, in a way similar to what happens to the Taylor shock profile solution of Burgers 
equation p. This, however, is still an open question. 
(iii)When Ac^/a^dpK^ = 1, then 

2c 
ad 

and the frequency is again a pure imaginary number, u = iuj, with 

This is a transition case between the two cases previously described, in which the system 
will again present either a strictly dissipative or anti-dissipative behavior, depending on 
whether uj is negative or positive. For these values of n, as well as for the values of case (i), 
the nonlinear dissipative Boussinesq system will present no oscillatory motion, which means 
that its solutions must describe pure nonlinear diffusive phenomena at this region. 



V. FINAL COMMENTS 

The results obtained in this paper are two- fold. First, from a perturbative analysis, 
corresponding to the long-wave in shallow-water approximation, we obtained a dissipative 
generalization of the Boussinesq system of equations as governing the bulk velocity and 
surface perturbations of a Benard-Marangoni phenomenon. This system of coupled evolution 
equations includes nonlinearity, dispersion and dissipation. The predominance of any one of 
them depends on the relation between e and 6. Three cases are of special interest. First, when 
M = —12, the dissipative term vanishes, and assuming that 5^ ~ e, the usual Boussinesq 
system of equations is obtained. Second, when the Marangoni number is far enough from the 
critical value, that is when M + 12 = 0{1), and assuming now that 5 ~ e, the dissipative 
term dominates over the dispersive one. Then, by neglecting the dispersive term we got 
Eq. (21), which is a strictly dissipative version of the Boussinesq system of equations. And 
finally, as an intermediate case, we may also consider the situation in which M + 12 = 0{6). 
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Assuming again that 5^ e, it results that the dispersive and dissipative terms of Eq. (15) 
are of the same order. In this case, no term is neglected, and the generalized Boussinesq 
system, Eq. (15), will govern the bulk velocity and surface perturbations of the Benard- 
Marangoni system. Upon specialization to waves moving, say, to the right, each one of these 
three cases will lead to a single evolution equation for the surface displacement, which will 
be, respectively, the KdV, Burgers and KdV-Burgers equations. It should be remarked that 
the main feature of the perturbative scheme adopted here was to allow for a clear distinction 
of the relative importance of each term in Eq.(15). By considering different relations between 
the parameters 6 and e, a criteria to collect terms of the dominant order, and to neglect 
higher order ones, was thus established. 

The second result of this paper refers to the appearance of a new equation, which we have 
called the strictly dissipative Boussinesq system of equations. It seems to be a nonintegrable 
equation, and the existence of analytical solutions is still an open question. Despite of this 
problem, we have performed here a linear analysis to get some insight on the nature of the 
physical phenomena it might describe. Finally, due to the fact that it shows up in a physical 
system, we believe this system of equation deserves further investigations as only few results 
exist with respect to dissipative systems. 
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APPENDIX: 

We give here a summary of the results described, but not explicitly shown in the text. 
The expressions for each term of the velocity expansion, already written in terms of uq and 
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Po, and considering only those terms that will contribute to the evolution equations up to 
orders e and 5^, are 



SR f 1 S \ 

S^R f 2 R R 

ua — -;rr ( — Poxxx + -tUqu H — rPoxt — ^iJ-oxxt 



24 



e S 

6^R' 



e5' 



1 



720 \ ~Poxtt 



The next terms of the expansion will not contribute to the evolution equations. In the same 
way, the only term of the pressure expansion that will contribute is 

P2 = -—POxx- 

From Eq. (3), with the corresponding boundary conditions, we can obtain expressions for 
9i, $3 and 9 5 in terms of Uq and po only: 

Ra ^I^a . . oR^o' 

Oi = -ed—uox + ed ^^(5(7 - ^)uoxt - d -^Poxx , 
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-Uoxt + ^5 



.Ra 



05 = e5^ 



/ 1\ 1 ,Ra .2 

1 + - Uoxt H POxx - O—Uoxtt - O—Uox 

\ aj ea 2 Ra 



Again, these are the only components of the temperature expansion that will contribute, up 
to orders e and 5^, to the evolution equations. 

Now, from Eq. (6), we can get po in terms of Uq and rj only: 



Po = e?7 - ^S—uox - e5^ 
K 



UOxt + [^ + ^]Vxx 



1 _ i? ■ 

"^UQxxx '^^Oxtt 



This equation allows us to rewrite the u^s, p's and ^'s in terms of uq and rj. Consequently, 
the system of evolution equations, Eq. (12), obtained from Eqs. (5) and (7) could also be 
written in terms of Uq and r] only. 
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